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Abstract 

In the framework of renormalon consideration in the approximation of a large 
number of quark flavours, a role of the anomalous dependence of gluon propagator 
on the scale p, is shown. 

In the QCD perturbation theory, an account for the next-to-leading order term over 
a s is connected with the uncertainty, caused by a choice of the energy scale, determining 
the aq(p 2 ) value, since the transition to the p scale leads to the substitution 

a s 0 2 ) « aq(p 2 ) (7 - a a (fi 2 )^ In , 

where /3q = (11 N — 2n/)/3, N is the number of colours, rif is the number of quark 
flavours. So, the physical quantity, represented in the form 

r = r 0 aq(p 2 ){ 1 + rqaq + 0(a 2 )} 

with the given order of accuracy, will get the form 

r = r 0 aq(p 2 )jl + (r x - ^-ln^ja s + 0(a 2 s )^ 

= r 0 a s (ji 2 ){l+fia s + O(a 2 s )}. 

Thus, the value of rq coefficient for the a s correction depends on the scale of aq(p 2 ) 
determination. In ref.[[lj, one offered to fix the p scale so that the rq coefficient does 
not contain terms, proportional to (3q. The latter procedure can be understood as one, 
leading to rq to be independent of the number of quark flavours rif ~ (3q. As was shown 
in refs.[@. such choice of scale has the strict sense in the framework of the 1 /rif 
expansion, where the a s correction to the gluon propagator is determined by the fermion 
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loop contribution into the vacuum polarization. So, the expression for such contribution 
depends on the regularization scheme, and in the next-to-leading order one has 


a s (g 2 )D(k 2 } n) 


a s (/i 2 ) 

k 2 



a s n f 

Qn 



(i) 


where D(k 2 ,n) determines the transversal part of the gluon propagator D^ b x (k 2 , //) = 
5 ab D(k 2 , /j>)(—g u \ + k v k\/k 2 ). The constant value C is defined by the renormalization 
scheme, so Cj^g = —5/3 and Cy = 0 in the so-called V-scheme jj]. 

The summing of (n/a s ) n contributions into the gluon propagator leads to the expres¬ 
sion 

a,^)D(k\^ = a ‘ ( // 2) , ( 2 ) 

i.e., in fact, it results in the account for the ’’r unning ” a s value in respect to the gluon 
virtuality. Note, the 1 /rif consideration is exact for the abelian theory, where a s D is 
the renormalization group invariant. Moreover, the e -c A qq D value does not depend on 
the renormalization scheme, and therefore it results in the scheme-independence for the 
expression 

j PoHePV/AhcD) ' 

As one can see in the performed consideration, the transition to the nonabelian 
theory is done by the 2rif/3 —> —(3q substitution, that is called as the procedure of 
’’naive nonabelization”, giving the correct ’’running” of the coupling constant. 

Further, consider the physical quantity r, for which the first order a s -contribution is 
calculated as an integral over the gluon virtuality with the weight F(k 2 ) 


r = / ^ ®s(g 2 )F(k 2 ) . 

Then the account for the ( rifa s ) n corrections to the gluon propagator leads to the sub¬ 
stitution of ’’running” value a s (e c k 2 ) for a s (g 2 ) 
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a s (k 2 )^o 
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( 3 ) 


so that in the offered procedure of the scale fixing [ij], the introduction of the first 
order n/a s -correction results in the substitution /i 2 —> ft 2 = /i 2 exp{C + (In k 2 /g 2 )}, 
where the average value is determined by the integral with the F weight. However, in 
some cases [{2|, |j] the integration of the n-th order term in expansion (||) leads to the n\ 
factorial growth of the coefficients for the expansion over a(?(/i 2 ), so that this rising can 
be determined by the region of low virtualities as well as large ones of gluon. Then one 
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says about the infrared and ultraviolet renormalons, respectively |2|. [5J, [4j. As was shown 
||], the renormalon leads to power-like uncertainties of the r evaluation, so 

(4) 

where k is positive for the infrared renormalon and it is negative for the ultraviolet 
one. For the two-point correlator of heavy quark vector currents, for instance, one has 
k — 4 and the corresponding uncertainty can be, in fact, eliminated in the procedure 
of definition for the nonperturbative gluon condensate, having the same power over the 
infrared parameter J3|. However, for the renormalized mass of heavy quark the infrared 
renormalon with k — 1 does not correspond to some definite quark-gluon condensate, and 
the value Am(/x) ~ A qcd can not be adopted into a definition of a physical condensate. 

In the present paper we modify the renormalon due to the account for the anomalous 
dependence of gluon propagator on the /i scale, and the modification turns out, basically, 
to be essential for numerical estimates, related with the determination of physical scale, 
fixing a s . 

In the covariant gauge, the gluon propagator has the form 


k 2) (- 9 ^ + “ «z(h 2 M^ 2 ))^) 


fc 2 u;(/i 2 , k 2 ) 


where aiu is the renormalization group invariant. The solution of one-loop equation of 
the renormalization group for a;(/i 2 ) in the MS -scheme allows one to represent the gluon 
propagator in the form || 


D? x (k?,n) = 


sab 


k 2 




s (m 2 )V 

s(Mq) ) 


~9v\ + (1 — b) 


Kkx 
k 2 


(5) 


where n a = (13N — 4n/)/(6/?o) and b is the arbitrary gauge parameter^. In eq.(^[) we take 
into account the anomalous dependence for the gluon propagator on the scale, u = a;(/i 2 ), 
only. 

The choice of fixed point uj = 1 corresponds to the standard renormalon. The 
consideration of the gluon propagator at u ^ 1 leads to the modified renormalon. 

In the leading order over 1 /rif one has n a = 1, and the cases with 1) u;(/ig) > 1, 
2) u;(/io) = 1 and 3) u(fi q) < 1 at some large /i 0 values can be formally come to the 
different choices of the renormalization group invariant /i 9 , such that uj(l^ 2 g) = 0 and 
1) n 9 < A qcd, 2) ji g = A qcd and 3) \i g > A QC d, respectively. The a s (/i 2 ) value is 
considered as the formal one-loop expression 
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A QCD 


1 In ref.p the solution with the fixed 6 = 3 value is considered and lo is expressed through 6/a;. 
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Thus, the 1 jnf consideration gives 


Qs(/j 2 ) _ Qs(yK(/ig) 

w(/i 2 ) a a (//2) - a,(/i 2 ) 


( 6 ) 


where d is the scheme-independent invariant of the one-loop renormalization group, and 
it is equal to 

j _ i ^-QCD _ 4 ?r 

“ n A 2 A)«,(a 2 ) ■ 

Further, one can repeat the calculation of (rifa s ) n contributions into the gluon propaga¬ 
tor with the substitution of value (^) for a s (p 2 ), so that this replacement corresponds to 
the account for the anomalous dependence of the gluon propagator versus the /i scale. 
Then the scale fixing due to the next-to-leading order a s correction results in the ex¬ 
pression 

jj 2 = p 2 exp{C + (In k 2 //i 2 ) + d} , 

and the summing of the corresponding contributions modifies eq.(|2|) 


a s (/r 2 )D(fc 2 ,/i) 


a s (e c+d k 2 ) 

V 2 


a s (e c k 2 ) 

uj(e c k 2 )k 2 


(7) 


Next, one can define the V -scheme, introducing Cy = —d. Then 


\V _ 5/W/2*AfS 

1V QCD ~ e 1V QCD ■ 

In the V -scheme the perturbative potential between the heavy quark and antiquark in 
the colour-singlet state will have the form 

= 4 (8) 

3 q- 

at q 2 = k 2 . Potential (||) comes to the Richardson potential ||, when one uses (q 2 + 
A q CD ) instead of a^(q 2 ), and A q C d fixes the linearly rising part of potential, confining 
quarks with the distance increase, 


AVUx) 


27 
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The fitting of mass spectra for the charmonium and bottomonium in the Richardson 
potential gives 


A^ cd = 400 ± 15 MeV . 
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The one-loop expression with a^' fS (rri 2 z ) = 0.117 ± 0.005 0 corresponds to A q^ d — 
85 ± 25 MeV, that allows one phenomenologically to determine the values 

e ~ d/ 2 = 2.1 ± 0 . 5 , 

nfs = i 8 o ±45 MeV, (9) 

a s (n 2 g ) = 0.94 ±0.21. 


However, an additional uncertainty into the estimates is involved by the transition from 
the number of flavours nf(mz ) = 5 to another value n/(l GeV) = 3. Obtaining eq.(§), 
one has supposed, that A qcd does not depend on the number of flavours, so that at 
the scale of ’’switching on (off)” the additional flavour of quarks = m nf+ i, the QCD 
coupling has a discontinuity, related with the step-like change /do(^/) — 1> Po( n f ± 1) = 
/ 3 0 (rif ) —2/3. One can avoid such discontinuities, if one supposes that A q C d depends on 
the flavour number, so that a s (fi 2 = +1 , Aq^ d , n f) = ay(/i 2 = m 2 /+1 , Aqq^\ rif + 1). 

Then in the one-loop approximation one finds 


2 

A (n f ) _ 

A QCD ~ A QCD 1 („ /+ i) ) 

a qcd 


( 10 ) 


Setting A.QQD = 85 ± 25 MeV, one gets 

e ~ d/ 2 = 1>24 ± 0>32 _ 


( 11 ) 


Hence, the a s rescaling, accounting for eq.(p±|), results in the d value, that agrees with 
the standard renormalon (d = 0) within the current accuracy. 

Next, one can use the nonabelian ’’running” of u with n a ^ 1 and the same scale 
as in the ’’running” of a s . In the latter procedure one has two scales Aqcd and /i 9 , so 
that two ’’measurements” of a s /u> quantity, entering physical values, are needed. If one 
fixes [i'y from the linear part of the heavy quark potential and uses the Z pole data for 
the Aqcd evaluation, then the obtained values show the standard renormalon validity 
within errors. 

Thus, we have found that in the 1 /rif consideration the account for the anomalous 
dependence of the gluon propagator on the fi scale results in the modification of renor¬ 
malon because of the presence of /i s scale in addition to A qcd, and the difference from 
the standard renormalon turns out to be essential in numerical estimates. We have also 
determined in a phenomenological manner. 
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